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This  paper  proposes  a  multi-choice  goal  programming  model  for  dealing  with  the  capacity  expansion 
planning  problem  of  the  renewable  energy  industry.  This  model  involves  decisions  regarding  the  optimal 
mix  of  different  plant  types,  location  selection  and  other  criteria.  Different  types  of  plants  should  be 
located  in  appropriate  places  so  as  to  minimize  the  total  deviations  from  predefined  goals  concerning 
power  generated,  investment  cost,  emission  avoided,  jobs  created,  operation  and  maintenance  costs, 
distance  security,  and  social  acceptance.  The  proposed  method  is  superior  to  the  goal  programming 
model  proposed  by  Ramon  and  Cristobal,  in  that  it  can  avoid  underestimation  of  aspiration  level,  expand 
the  potential  feasible  region,  and  achieve  findings  more  closely  approach  actual  conditions.  In  addition, 
the  social  acceptance  of  the  renewable  energy  planning  problem  in  Taiwan  is  modeled  by  the  MCGP  to 
demonstrate  its  usefulness. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Because  of  continuous  economic  growth  and  industrialization 
worldwide  and  the  soaring  energy  demands  in  the  five  major 
developing  countries  of  Asia  (China,  India,  Indonesia,  Korea  and 
Thailand),  the  demands  on  fossil  energy  for  industry,  transporta¬ 
tion  and  daily  life  are  becoming  increasingly  large.  Table  1  shows 
the  global  reserves  and  the  projected  number  of  years  of  supply 
remaining  for  major  fossil  fuels  [1,2].  Even  with  recent  discoveries 
of  new  reserves,  such  as  shale  gas,  tight  sands  and  shale  oil  in 
North  America,  as  well  as  shale  oil  in  many  other  areas,  fossil  fuels 
are  gradually  being  depleted.  International  oil  and  coal  prices 
exceeded  US  $102  per  barrel  and  US  $75  per  ton  in  April  2010  [3,4]. 
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Since  Taiwan  has  almost  no  energy  resources,  97.92%  of  the  total 
energy  requirement  depended  on  imports  in  2013  and  almost  all  of 
the  fossil  fuel  is  imported  from  politically  unstable  areas,  such  as 
the  Middle  East,  so  Taiwan's  energy  security  has  become  increas¬ 
ingly  vulnerable.  Energy  imports  also  increased  from  3.88%  of 
Taiwan's  GDP  in  2002  to  14.55%  in  2012,  which  creates  a  large 
drain  on  Taiwan's  economy.  In  addition,  traditional  uses  of  fossil 
fuels  aggravate  environmental  pollution  and  the  greenhouse  effect. 
Thus,  the  development  of  renewable  energy  is  an  increasingly 
important  issue  for  both  developing  and  developed  countries. 

Taiwan's  energy  supply  and  consumption  in  2012  are  listed  in 
Fig.  1.  Compared  with  2011,  petroleum  is  3.4%  higher,  renewables 
are  0.1%  higher,  nuclear  is  0.2%  higher,  natural  gas  is  0.7%  higher, 
and  coal  is  4%  lower.  This  means  that  Taiwan  is  gradually  improving 
with  the  use  of  cleaner  fuels  to  reduce  carbon  emission  from  power 
generation.  Since  the  2011  nuclear  disaster  at  Fukushima,  the 
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Table  1 

Global  reserves  and  availability  of  major  energy  resources. 


Category  item 

Oil  (  +  oil  sanss) 

Natural  gas 

Coal 

Total  reserves  (end  of  2011) 

1481  billion  barrels 

196  trillion  cubic  meters 

860  billion  tons 

Yield  (2011) 

30  billion  barrels 

3.2  trillion  cubic  meters 

6.7  billion  tons 

Available  years 

49  years 

61  years 

112  years 

Fig.  1.  2012  Energy  supply  and  consumption  in  Taiwan. 

government  of  Taiwan  has  considered  a  nuclear  energy  reduction 
policy.  Overall,  it  is  a  top  priority  of  Taiwan  to  develop  clean, 
sustainable,  and  independent  energy  and  achieve  a  balance  between 
energy  security,  environmental  protection,  industrial  competitive¬ 
ness,  and  economic  development.  This  is  particularly  important  for 
Taiwan's  environment,  with  the  sunny  weather  in  southern  Taiwan 
and  abundant  wind  resources  on  the  west  coast  of  Taiwan.  These 
issues  lead  to  both  crisis  and  opportunity  for  Taiwan  in  becoming  a 
sustainable  community.  In  order  to  encourage  the  development  of 
renewable  energy,  the  National  Science  Committee  (NSC)  initiated 
the  flagship  project  “New  Projects  of  Energy  Technology”  in  2008. 
Subsequently,  the  Ministry  of  Economic  Affairs  (MOEA)  promulgated 
the  “Renewable  Energy  Development  Act  (REDA)”  in  2009,  the 
“Incentive  Program  of  Offshore  Wind  Power  Demonstration  System” 
in  2012,  the  “Million  Solar  Rooftop  Photovoltaic  Program”  in  2013, 
and  the  “Small  Wind  Turbine  Generation  System  Demonstration 
Incentives”  in  2014.  Many  kinds  of  financial  compensation  are  offered 
by  the  central  and  local  governments.  Further  programs  are  also 
expected  under  the  Bureau  of  Energy  (BOE)  for  promoting  renewable 
energy  usage:  (i)  Total  renewable  energy  capacity  of  9952  MW  in 
2025  (i.e.,  15.5%  of  the  total  power  generation  installation  capacity  in 
Taiwan);  (ii)  1000  on-  and  off-shore  wind  turbines,  with  the  wind 
power  capacity  of  4200  MW  in  2030  and  progress  on  the  construc¬ 
tion  of  the  wind  power  system  as  shown  in  Table  2;  and  (iii)  the 
“Million  Solar  Rooftop  Photovoltaic  Program”  project,  with  photo¬ 
voltaic  power  capacity  of  3100  MW  in  2030.  In  particular,  the 
Industrial  Technology  Research  Institute  (ITRI),  Taiwan  Power  Com¬ 
pany,  Chinese  Petroleum  Corporation,  and  others  have  begun  to  carry 
out  research  on  renewable  energy  technology.  From  the  promulga¬ 
tion  of  the  REDA  in  2009  to  the  end  of  2012,  the  capacity  of 
photovoltaic  power  increased  from  11  MW  to  222.4  MW,  the  capa¬ 
city  of  wind  power  increased  from  436  MW  to  621  MW. 


However,  these  infrastructure  projects  are  often  contested  and 
face  public  resistance,  in  particular  against  many  of  the  on-shore 
wind  power  projects  in  Taiwan.  Therefore,  public  resistance  should 
be  considered  as  an  important  factor  in  a  model  for  the  optimal 
location  of  renewable  energy  facilities.  The  capacity  expansion 
planning  problem  of  the  renewable  energy  industry  should  not 
only  to  determine  the  best  renewable  energy  project  among 
different  alternatives,  but  also  develop  expansion  plans  for  energy 
generation,  transmission  and  distribution  systems  with  considera¬ 
tion  of  the  potential  public  resistance  to  these  projects.  This 
problem  involves  important  decisions  regarding  the  optimal  solu¬ 
tions  of  different  plant  types,  locations,  and  capacities.  A  variety  of 
approaches  have  been  employed  to  model  the  generation  and 
distribution  of  electricity,  such  as  of  multi-criteria  decision-mak¬ 
ing  models  [5-7],  optimization  models  with  geographic  informa¬ 
tion  systems  for  regional  renewable  energy  development  [8],  and 
goal  programming  (GP)  model  for  the  optimal  mix  and  location  of 
renewable  energy  plants  [9  .  As  a  result,  multi-criteria  decision¬ 
making  methods  such  as  the  Analytical  Hierarchy  Process  (AHP), 
Elimination  and  Choice  Expressing  Reality  (ELECTRE),  Preference 
Ranking  Organization  Method  for  Enrichment  of  Evaluations 
(PROMETHEE),  Technique  for  Order  Preference  by  Similarity  to 
an  Ideal  Solution  (TOPSIS),  and  VIsekriterijumsko  KOmpromisno 
Rangiranje  (VIKOR)  have  been  extensively  adopted  for  the  selec¬ 
tion  of  renewable  energy  problems  such  as  wind-farm  projects, 
geothermal  projects  and  hydro-site  selection  [10-16]. 

Though  these  methods  have  been  widely  applied  to  the  renew¬ 
able  energy  industry,  the  multi-choice  goal  programming  (MCGP) 
method  [17-19]  has  rarely  been  used  for  the  location  selection 
problem  of  renewable  energy  plants  and  capacity  expression 
decisions.  MCGP  is  a  multi-choice  aspiration-level  approach  which 
has  been  applied  to  many  real-world  problems,  such  as  housing 
selection  [20],  low-cost  carriers'  networks  [21],  supplier  selection 
[22,23],  forestry  plans  [24],  forestry  management  [25]  and  a 
piloted  quality-management  system  26].  Moreover,  using  MCGP 
to  deal  with  the  optimal  location  of  renewable  energy  plants,  can 
provide  the  following  advantages: 

•  MCGP  expands  the  original  feasible  region  (it  may  be  obtained 
using  Weight  GP,  MINMAX  GP,  or  Lexicographic  GP)  to  the 
potential  feasible  region  for  obtaining  the  appropriate  solution 
for  decision-makers  (DMs)  responsible  for  the  selection  of  the 
renewable  energy  plant  location  17]. 

•  It  not  only  can  deal  with  quantitative  issues  (e.g.,  minimization 
of  operation  and  maintenance  costs)  but  also  qualitative 
problems  (e.g.,  maximization  of  user  satisfaction  and  social 
acceptance) [17  . 

•  It  allows  DMs  to  set  multiple  aspiration  levels  to  avoid  under¬ 
estimating  settings  determined  by  conservative  chief-executive 
officers.  For  example,  the  lower  power  capacity  of 
113,609,530  W  is  obtained  in  the  model  of  Ramon  and  Cristobal 
(RC)  [9]  because  the  power  generation  goal  (110  x  106  W)  is 
underestimated  in  their  model.  In  contrast,  a  higher  power 
capacity  of  116,587,600  W  is  obtained  by  the  MCGP  model 
because  MCGP  can  reflect  higher  aspirations.  That  is,  by 
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avoiding  underestimation  of  the  aspiration  level,  it  is  possible 
to  have  better  solution  in  MCGP  [19]. 

This  study  uses  MCGP  to  deal  with  the  selection  of  renewable 
energy  plants,  locations  and  the  capacity  expression  problem  by 
taking  into  account  the  above  advantages  and  striving  to  achieve 
more  reality-based  findings.  The  superiority  of  the  proposed 
method  can  be  seen  in  the  comparison  of  the  solutions  between 
Ramon  and  Cristobal  [9]  and  the  proposed  MCGP  method.  In 
addition,  the  MCGP  is  used  to  solve  a  wind  power  plant  location 
selection  problem  in  Taiwan  taking  into  account  public  resistance. 

2.  Ramon  and  Cristobal  model  (RC  model) 

In  order  to  locate  five  renewable  energy  plants  for  power 
generation  in  five  places  located  in  the  region  of  Cantabria  in  the 
north  of  Spain,  Ramon  and  Cristobal  [9]  proposed  a  novel  GP 


model  to  address  the  problem  of  how  to  locate  seven  candidate 
plants  (alternatives)  in  seven  locations:  La  Braguia,  La  Vega, 
Estacas,  Lunada,  Potes,  Reinosa  and  Santander.  The  attributes 
considered  when  evaluating  these  renewable  energy  systems  in 
this  model  are:  power  generated  (W),  investment  cost  (I),  tons  of 
emissions  of  C02  avoided  per  year  (TC02/y),  jobs  created  (J), 
operation  and  maintenance  costs  (OM),  distance  between  plants 
(D)  and  social  acceptance  (S).  Table  3  shows  the  data  for  the  first 
five  attributes.  Table  4  shows  the  distance  between  the  locations. 
Table  5  shows  the  social  acceptance  by  the  local  population  for  the 
hypothesized  realization  of  the  renewable  energy  projects.  Social 
acceptability  is  expressed  using  a  scale  of  1  (low  acceptance)  to  10 
(high  acceptance).  As  different  types  of  plants  can  be  placed  in 
each  location,  the  goal  is  to  ensure  the  minimum  total  deviation 
from  the  goals.  These  goals  are  given  as 

-  (Gl)  The  power  generated  must  be  higher  than  110  x  106  W. 

-  (G2)  The  investment  cost  must  be  limited  to  100  x  106€/year. 


Table  2 

1000  wind  turbines  project  in  Taiwan. 


Years 

2012 

2015 

2020 

2025 

2030 

On-shore  wind 

621 

866 

1200 

1200 

1200 

(MW) 

(314  turbines) 

(350  turbines) 

(450  turbines) 

(450  turbines) 

(450  turbines) 

Off-shore  wind 

0 

15 

600 

1800 

3000 

(MW) 

(4  turbines) 

(120  turbines) 

(360  turbines) 

(600  turbines) 

Total 

621 

881 

1800 

3000 

4200 

(314  turbines) 

(354  turbines) 

(570  turbines) 

(810  turbines) 

(1025  turbines) 

Table  3 

Alternatives  for  electric  generation. 

Alternative 

W 

I 

TC02/y 

J 

OM 

Wind  power  (10  <  P  <  50  MW) 

60,959,000 

23,425,000 

9,649,680 

15 

37,750 

Hydroelectric  (P  <  10  MW) 

3,940,100 

7,500,000 

472,812 

8 

7250 

Hydroelectric  (10  <  P  <  25  MW) 

962,000 

14,000,000 

255,490 

8 

14,000 

Hydroelectric  (25  <  P  <  50  MW) 

412,000 

21,035,000 

255,490 

12 

21,000 

Solar  (P>  10  MW) 

306,031 

5,320,000 

482,856 

10 

1792 

Biomass  (P  <  5  MW) 

13,612,500 

9,015,000 

2,524,643 

15 

27,100 

Biomass  (P  >  50  MW) 

37,770,000 

47,958,400 

4,839,548 

20 

75,000 

Table  4 

Distance  between  locations. 

Alternative 

La  Braguia 

La  Vega 

Estacas 

Lunada 

Potes 

Reinosa 

Santander 

La  Braguia 

— 

6.5 

15 

39 

126 

61 

53 

La  Vega 

- 

8 

32 

121 

54 

65 

Estacas 

- 

24 

128 

62 

72 

Lunada 

- 

143 

78 

56 

Potes 

- 

130 

107 

Reinosa 

- 

75 

Santander 

— 

Table  5 

Social  acceptance  (1-10)  of  each  plant  in 

each  location. 

Alternative 

La  Braguia  (8) 

La  Vega  (9) 

Estacas  (10) 

Lunada  (11) 

Potes  (12) 

Reinosa  (13) 

Santander  (14) 

( 1 )  Wind  power  ( 10  <  P  <  50  MW) 

7 

9 

7 

9 

(2)  Hydroelectric  (P  <  10  MW) 

5 

6 

5 

7 

(3)  Hydroelectric  (10  <  P  <  25  MW) 

6 

7 

6 

7 

(4)  Hydroelectric  (25  <  P  <  50  MW) 

(5)  Solar  (P>  10  MW) 

4 

5 

8 

7 

9 

6 

(6)  Biomass  (P<5  MW) 

9 

8 

8 

9 

(7)  Biomass  (P>50  MW) 

9 

8 

9 
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-  (G3)  The  emissions  avoided  must  be  higher  than  18  x  106  t. 

-  (G4)  The  jobs  created  must  be  higher  than  70. 

-  (G5)  The  operation  and  maintenance  costs  must  be  limited  to 
350,000€/year. 

-  (G6)  The  distance  between  plants  must  be  maximized  to 
2910  Km. 

-  (G7)  The  social  acceptance  must  be  as  close  as  possible  to  the 
highest  level  of  50. 

Fig.  2  shows  the  multi-source  multi-sink  network  involving 
seven  alternatives  and  seven  locations,  a  flow  of  k  units  from  the 
source  node  through  alternatives  and  place  nodes  to  the  sink  node 
represents  an  assignment  of  k  plants  to  places  in  which  k 
compatible  couples  are  created.  This  behavior  can  ensure  that 
alternatives  and  place  nodes  are  represented  by  binary  variables. 
Then  Ramon  and  Cristobal  [9]  proposed  a  GP  model  (see  [9]  for 
details)  for  the  optimal  mix  and  location  of  renewable  energy 
plants.  In  their  model,  the  optimal  mix  of  different  plant  types, 
locations,  and  capacity  expansion  can  be  easily  determined  using 
the  multi-source  multi-link  network  analysis  mechanism.  The 
compromise  solutions  of  RC  model  thus  obtained  are  shown  in 
Table  6. 

3.  Improvement  of  application 

The  GP  approach  is  a  well-known  and  popular  multi-objective 
mathematical  programming  model.  It  can  be  accomplished  with 
various  types  of  methods  such  as  Lexicographic  GP,  Weight  GP, 
MINMAX  (Chebyshev)  GP,  and  Extended  GP  [27  .  Its  inherent 
flexibility  allows  users  to  take  into  account  several  conflicting 
objectives,  multiple  criteria,  and  incomplete  information  at  same 

Alternatives  Places 


time  for  choosing  the  most  satisfactory  solution  within  a  feasible 
region.  GP  was  first  introduced  by  Charnes  et  al.  (1955)  and  further 
developed  by  Charnes  and  Cooper  (1961),  Lee  (1972),  lgnizio 
(1985),  Tamiz  et  al.  (1998)  and  Romero  (2001),  among  others 
[17].  In  practice,  the  GP  model  minimizes  2,(dj+  +df ),  subject  to 
constraints  /j(x)-dj+  -pdf  =gh  where  xeXcRn,  d,+  and  df  are 
positive  and  negative  deviations  of  the  achievement  function  /,(x) 
from  the  scalar  aspiration  level  g{.  In  other  words,  the  purpose  of 
GP  is  to  minimize  the  deviations  between  the  achievement  of 
goals  and  the  scalar  aspiration  levels  (it  may  be  predefined  by  DMs 
for  their  problems). 

In  contrast,  the  role  of  MCGP  is  to  aid  a  DM  in  expanding  the 
original  feasible  region,  which  may  be  obtained  by  GP,  to  the 
potential  feasible  region  to  find  more  and  better  solutions. 
The  MCGP  model  minimizes  ,(df+ +df +ef +ef ),  subject  to 
constraints  /,(*) -  d,+  +d~  =  yh  yt - e+  +  e ~  =  g,  max  or  g,  min , 
where  x  eX  c  Rn,  df+  and  d,  are  positive  and  negative  deviations 
of  the  achievement  function  /,(x)  from  the  vector  aspiration  level 
yb  ef  and  er  are  positive  and  negative  deviations  of  the  vector 
achievement  level  y{  from  the  scalar  aspiration  level  g,  max  or 
gu  min.  Here,  gu  max(g,  min)  is  the  sufficiently  large  value  for  the  case 
of  the  more  (less)  the  better  in  real  problems.  In  order  to  more 
closely  approach  reality  with  a  potential  solution  than  the  com¬ 
promise  solution  obtained  by  Ramon  and  Cristobal  [9  ,  the 
aspiration  for  each  goal  in  the  model  of  Ramon  and  Cristobal 
should  be  changed  as  follows: 

-  (Gl)  Power  generated  must  be  higher  than  110  x  106  W  and 
the  higher  the  better. 

-  (G2)  Investment  cost  must  be  limited  to  100  x  106€/year  and 
the  lower  the  better. 

-  (G3)  Emissions  avoided  must  be  higher  than  18  x  106  t  and  the 
higher  the  better. 

-  (G4)  Jobs  created  must  be  higher  than  70  and  the  higher  the 
better. 

-  (G5)  Operation  and  maintenance  costs  must  be  limited  to 
350,000€/year  and  the  lower  the  better. 

-  (G6)  Distance  between  plants  must  be  maximized  to  2910  Km 
and  the  greater  the  better. 

-  (G7)  Social  acceptance  must  be  as  close  as  possible  to  the 
highest  level  of  50  and  the  higher  the  better. 


Table  7 

Places  to  locate  the  plants  of  various  models. 


Alternative 

MCGP-MSMS  model 

( 1 )  Wind  power  ( 10  <  P  <  50  MW) 

Estacas  (10) 

(2)  Hydroelectric  (P  <  10  MW) 

Reinosa  (13) 

(5)  Solar  (P>  10  MW) 

Potes  (12) 

(6)  Biomass  (P<5  MW) 

Lunada  (11) 

(7)  Biomass  (P>50  MW) 

Santander  (14) 

Table  6 

Solutions  comparison  of  various  models. 


W  (Gl) 

I  (G2) 

TC02/y  (G3) 

J(G4) 

OM  (G5) 

Distance  (G6) 

Social  acceptance  (G7) 

Goals 

110  x 106 

100  x 106 

18  x  106 

70 

350,000 

2910 

50 

RC  model 

113,609,530 

99,718,400 

17,752,217 

68 

155,642 

2323 

43 

MCGP-MSMS  mdel 

116,587,600 

93,218,400 

17,969,540 

68 

148,892 

1043 

42 

MCGP-MSMS-normalized 

116,587,600 

93,218,400 

17,969,540 

68 

148,892 

1043 

42 

MCGP-MSMS-normalizeda 

113,059,500 

106,753,340 

17,752,220 

72 

162,642 

1043 

41 

MCGP-MSMS-normalized 

116,587,600 

93,218,400 

17,969,540 

68 

148,892 

948 

43 

a  10  times  as  a  weight  add  to  job  created  goal. 
b  10  times  as  a  weight  add  to  social  acceptance  goal. 
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For  treating  these  revised  goals  with  new  aspiration  levels,  this 
study  proposes  a  novel  method,  referred  to  as  the  MCGP  with  a 
multi-source  and  multi-sink  network  (MCGP-MSMS)  model  for 
solving  the  capacity  expansion  planning  problem  of  the  renewable 
energy  industry  (see  the  Appendix  for  details).  This  MCGP-MSMS 
model  is  then  solved  using  LINGO  [28]  to  obtain  the  solution,  as 
shown  in  Tables  6  and  7.  Table  6  shows  the  deviations  from  the 
goals,  while  Table  7  and  Fig.  3  show  the  selected  locations  for  each 
plant.  As  seen  in  Table  6,  goals  Gl,  G2  and  G5  are  achieved,  while 
goals  G3,  G4,  G6  and  G7  are  not  in  both  RC  and  MCGP-MSMS 
models. 

It  is  noteworthy  that  the  power  generated  (Gl),  investment 
cost  (G2),  emissions  avoided  (G3),  and  operation  and  maintenance 
cost  (G5)  achieved  in  MCGP-MSMS  model  are  better  than  the 
results  obtained  by  the  model  of  RC  [9].  Compared  with  the  RC 
model,  the  power  generated  is  2.62%  higher,  investment  cost  is 
6.97%  lower,  emissions  avoided  are  1.22%  lower,  and  operation  and 
maintenance  costs  are  4.53%  lower  in  the  MCGP-MSMS  model. 
Obviously,  the  underestimation  problem  can  be  overcome  by  the 
proposed  MCGP-MSMS  model  and  the  original  feasible  region  in 
the  GP  has  been  expended  to  the  larger  potential  feasible  region  in 
MCGP-MSMS  [17].  The  features  of  the  MCGP-MSMS  and  GP  are 
shown  in  Fig.  4. 


Alternatives  Places 


/. 


Deviational  variables  df  t  d[~ ,  ef+  and  er  in  the  model  of  RC  [9] 
or  MCGP-MSMS  can  vary  in  different  ranges,  resulting  in  an 
unintentional  bias  towards  objectives  with  a  large  magnitude. 
For  example,  the  goals  of  power  generated  and  investment  cost  in 
the  model  of  RC  [9]  require  a  relatively  large  value  to  obtain  a  large 
value  in  the  objective  function.  This  forces  these  two  goals  to  be 
given  higher  priority  than  other  goals  to  reach  the  target  in  the 
optimization  model.  Therefore,  the  use  of  a  normalization  techni¬ 
que  to  overcome  this  problem  is  an  important  issue  in  multiple- 
criteria  decision-making.  In  order  to  solve  this  problem,  the  utility 
normalization  technique  (Chang  17] )  is  slightly  modified  as 
follows: 

Let  nk  =  [0,Uk]  and  pk  =  [0,Vk]  Wk  where  Uk  and  Vk  are  the 
upper  bounds  of  nk  and  p]o  respectively.  Dk  =  Max  { Uk ,  Vkj  V  k  and 
Lj  =  Min  [Uk,  Vk }  Vk.  The  normalized  weights  are  easily  calcu¬ 
lated  as  wk  =  (Lj/Lj+Dk),  Wj  =  (Dk/Lj+Dk ),  and  wi  =  (wjUk/Ui) 
(i  e  {1,  ...,I<}\\j}).  This  ensures  that  wk  (k=  1,2,  ...K)  can  have 
roughly  the  same  magnitude  as  wkUk  =  WjU;  =  WjUj. 

For  example,  assume  that  the  deviational  variables  dk  and  nk  in 
the  MSCP-MSMS  model  can  vary  in  different  ranges  as 
d\,rii  =[0,120  x  10s],  d2,n2  =  [ 0,100  x  10s],  d3,n3=[ 0,19  x  10s], 
d4,  n4  =  [0, 90],  d5,n5  =  [0,350000],  d6,  n6  =  [0,4000],  and  d7,n7  = 
[0, 50],  resulting  in  unintentional  bias  toward  objectives.  According 
to  the  above  idea,  the  values  of  Uk  and  Vk  can  be  computed  as 
follows: 

Lb  =  V,  =  120  x  106,  U2  =  V2  =  100  x  106,  U3  =  V3  =  19  x  106, 
V4  =  V4  =  90,  U5=V_ 5  =  350000,  U6  =  V6  =  4000,  jind  U7  = 

V7  =  50.  Dk  =  max  {Uk,  Vk)  =  120  x  106,  Lj  =  min  {Uk,  Vk)  =  50, 
W\  =  (Lj/Lj  -\-Dk)  =  0.41 666  x  10  w7  =  (Dj</Lj-\-Dk)  =  0.999999, 
w2  =  (u7U7/U2)  =  0.499999  x  10~6,  w3  =  (u7U7/U3)  =  0.2631576 
xlO  W4  =  (u7U7/U4)  =  0.555555,  W5  =  (u7U7/U3)  =  0.142857 
xl0~3,  and  w6  =  ( u7U7/U6 )  =  0.0124999. 

This  ensures  that  dk  and  nk  (/<  =  1,2, ...,  7)  can  have  roughly  the 
same  magnitude  as  w^U-i  =  ...  =  w7U7  «  50  for  deviational  vari¬ 
ables  dk  and  nk  in  the  objective  function  of  the  MSCP-MSMS 
model.  The  normalized  weights  applied  to  the  MCGP-MSMS 
model  can  be  expressed  as  follows: 

(MCGP-MSMS-normalized) 

Min  (df  +ef )+w2(d2h  +e  +  )+w3(d3_  +e3-)+w4(d4  +e^)  + 

w5(d5+  +e+)+w6(d6“  +eg  )+w7df 

s.t.  Eqs.  (3.1)-(3.36). 

The  MCGP-MSMS-normalized  model  is  again  solved  using 
LINGO  [28]  to  obtain  the  same  solution  as  in  the  MCGP-MSMS 
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Fig.  4.  Possible  feasible  regions  in  various  GP  methods. 


Fig.  5.  The  solution  of  MCGP-MSMS-normalized  with  job  creation. 
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Alternatives 


Places 


Fig.  6.  The  solution  of  MCGP-MSMS-normalized  with  social  acceptance. 


model.  In  order  to  demonstrate  the  usefulness  of  the  proposed 
normalized  mechanism,  consider  the  following  cases: 

Case  I:  owing  to  the  increased  unemployment  rate  in  Spain, 
DMs  hope  to  increase  the  employment  rate,  so  the  weight  added 
to  the  goal  of  job  creation  is  increased  by  10  times  in  the  MCGP- 
MSMS-normalized  model.  Solving  this  problem  with  LINGO  [28] 
yields  a  new  solution,  as  shown  in  Table  6  and  Fig.  5.  As  can  be 
seen  in  Table  6,  Goals  Gl,  G4,  and  G5  are  achieved,  while  goals  G2, 
G3,  G6  and  G7  are  not  in  both  RC  and  MCGP-MSMS-normalized 
models.  Notably,  a  value  of  72  for  the  goal  of  job  creation  is 
achieved,  while  the  highest  achievements  of  other  goals  decline 
slightly.  It  is  clear  that  using  the  normalized  mechanism,  DMs  can 
define  the  appropriate  weights  for  the  MCGP-MSMS-normalized 
model  to  achieve  their  particular  goals. 

Case  II:  due  to  public  resistance,  DMs  would  like  to  increase  the 
social  acceptance  rate  to  the  maximum  possible  value,  and  to  do  so 
the  weight  added  to  the  goal  of  social  acceptance  is  increased  by  10 
times  in  the  MCGP-MSMS-normalized  model.  This  problem  is 
solved  with  LINGO  [28]  to  obtain  the  solution  shown  in  Table  6 
and  Fig.  6.  Fig.  6  shows  the  locations  for  each  plant.  As  seen  in 
Table  6,  Goals  Gl,  G2,  and  G5  are  achieved,  while  goals  G3,  G4,  G6 
and  G7  are  not  in  both  RC  and  MCGP-MSMS-normalized  models. 
Notably,  a  value  of  43  for  the  goal  of  social  acceptance  rate  is 
achieved,  which  is  the  highest  social  acceptance  rate  in  this  case, 
while  the  highest  achievements  of  other  goals  are  slightly  lower. 
Compared  with  the  RC  model,  the  power  generated  is  2.62%  higher, 
investment  cost  is  6.97%  lower,  emissions  avoided  are  1.22%  higher, 
and  operation  and  maintenance  costs  are  4.53%  lower  in  the  MCGP- 
MSMS-normalized  model.  It  is  seen  that  by  adding  weights  to  the 
objective  function,  the  proposed  method  can  easily  be  used  as  a 
decision  aid  to  determine  the  best  or  most  appropriate  solution  to 
multiple  objective  problems.  Each  goal  of  the  multiple  objective 
problem  can  be  divided  to  multiple  aspiration  levels  to  better  suit 
management  requirements,  such  as  “the  more  the  better  “  or  “the 
less  the  better”.  These  constraints  can  also  be  easily  added  in  the 
MCGP-MSMS-normalized  model  to  mirror  real-world  situations. 
The  MCGP-MSMS-normalized  model  provides  a  feasible  and  robust 
way  to  choose  an  optimal  location  for  renewable  energy  plants. 


4.  Social  acceptance  of  a  wind  turbine  installation 

Due  to  the  abundant  wind  resources  on  the  west  coast  of  Taiwan, 
28  possible  on-shore  locations  are  evaluated  for  implementation  of 


wind  turbine  generation  systems,  as  shown  in  Fig.  7.  Though  there 
were  314  wind  turbine  systems  installed  in  2012,  generating  621  MW, 
due  to  the  greater  awareness  of  environmental  protection  issues  and 
participatory  democracy,  protests  against  proposed  wind  turbine 
installation  have  occurred,  such  as  the  public  resistance  from  Yuan  Li 
township,  Penghu  county,  XinWu  township,  Tongxiao  township, 
Sanzhi  district,  Zhunan  township,  and  others.  As  shown  in  Table  2, 
there  are  886  wind  turbines  scheduled  to  be  installed  by  2030,  though 
due  to  popular  protests  this  will  be  the  most  difficult  part  of  the 
subsequent  installation  of  wind  turbines  in  Taiwan.  In  the  past  decade, 
there  have  been  many  studies  on  quantitative  goals  (e.g.,  maximize 
power  generated)  in  determining  the  optimal  location  of  renewable 
energy  facilities.  However,  qualitative  goals,  such  as  individual  pre¬ 
ferences  (e.g.,  minimize  popular  protests)  have  rarely  been  addressed 
in  solutions  to  this  problem.  The  second  objective  of  this  study  is  to 
provide  a  feasible  method  for  dealing  with  the  qualitative  problem  of 
reducing  public  resistance,  in  order  to  facilitate  subsequent  installation 
of  wind  power  system.  For  simplicity,  the  linear  utility  function  is 
given  as  the  following  two  membership  functions  (MF),  to 

represent  social  acceptance  rate  for  the  wind  turbines  installation.  A 
graph  representation  of  MF  is  shown  in  Fig.  8.  (4.1)  is  used  to 
maximize  /,(x)  >  g,  in  order  to  increase  the  social  acceptance  rate 
(e.g.,  the  greater  the  financial  compensation,  the  higher  the  social 
acceptance  rate). 


1, 

fj(x)  -  lj 


gi-li  ’ 

0, 


if  ffx )  >  gi 

if  k  <fi(x )  <  gi ,  (the  right  linear  MF) 
if  ffx)  <  k 


(4.1) 


(4.2)  is  used  to  minimize  /,•(*)  <gk  in  order  to  increase  the  social 
acceptance  rate  (e.g.,  the  less  the  wind  turbine  noise,  the  higher  the 
social  acceptance  rate). 


1. 


t/i</i(x»  =  < 


Uj  -//(*) 

-i 

0, 


Ui  ~gi 


V  ffx)  <  gi 

if  gi  <ffx)  <  u, ,  (the  left  linear  MF) 
if  ffx)  >  u, 


(4.2) 


where  u,  and  are  the  respective  upper  and  lower  tolerance  limits  for 
the  ith  goal.  The  achievement  model  with  the  MFs,  can  be  formulated 
as  follows: 


(MCGP-U) 

Min  ^  [cq(dj+  +  dj  )+/?j(e^+  -t-£j  )+2 /,•  -t ] 

i  =  1 


\ 

s.t.  rjfffx))  i/‘>  the  right  linear  MF 


ffx)-dl+  +d,. 
>Ji(ffX))+fi~ 


=  gh  the  upper  limit  value  g, 
=  1,  the  highest  value  of  1 


for  ffx)  >  gt 


/ 


rjj(fi(x))  =  ui ,  lix\  for  left  linear  MF 

11  1  Llj  gj 

fi(x )  -  ef+  +  er  =  gh  the  lower  limit  value  g, 
m(fi(x))  +  K  =  1>  the  highest  value  of  1 


\ 


for  ffx)  <  gk 


/ 


d,+ ,  dj  ,fi  ,  Hiiffx))  >  0 
x  e  F  (F  is  a  feasible  set) 

where  a,, ^,,2,,  and  y/i  are  the  weights  attached  to  the  deviational 
variables,  dt+  ,d[  ,ef+  ,ef  ,/f ,  and  hf . 

In  general,  Taiwan  citizens  are  opposed  to  high-polluting  coal- 
fired  and  fossil  fuel  power  generation,  and  the  safety  of  nuclear 
power  is  also  an  issue  of  public  concern,  so  there  is  support  for 
renewable  energy.  But  residents  still  do  not  want  to  have  wind 
turbines  sited  too  closed  to  their  dwellings,  schools  or  workplaces 
since  wind  turbines  are  a  source  of  noise  in  rural  environments. 
The  most  common  complaints  and  source  of  social  resistance  to 
wind  farms  are  noise  and  blade  flicker  (glare  reflecting  off  the 
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Kinmen  Jinsha 
Power  Plant  4MW*2 


Dayuan  Guanyin 
Power  Plant 
30MW*20 


Linkou  Power 
Plant  6MW*3 


Shihmen  Power 
Plant  3.96MW*6 


Tatan  Power 
Plant  15.1MW*8 


f - 

Siangshan  Power 
Plant  12MW*6 


Taichung  Port  Power 
Plant  34MW*17 


Taichung  Power 
Plant  8MW*4 


Changhua  Coastal 
Industrial  Park  Power 
Plant  62MW*31 


Wanggong  Power 
Plant23MW*10 


Mailiao  Power 
Plant  46MW*23 


Yunlin  Sihhu 
Power  Plant 
28MW*  14 


Penghu  Power 
Plant 

10.2MW*14 


Guanwei  Taowei  Power 
Plant43.7MW*19 


Chuenfong  Power 
Plant  1.75MW*1 


Sinfonq  Power 


PlantlT.5MW*5 


Chyidiing  Jwunan  Power 
Plant14.7MW*7 


Lungwei  Power 
Plant  34MW*16 


Dapeng  Power 
Plant  42MW*21 

Lunggang  Power 
Plant  11.5MW*5 

Taichung  Daan  Dajia 
Power  Plant 
75.9MW*33 


Luhwei  Power 
Plant 

96.6MW*42 


Formosa  Heavy 
Industries  Mailiao 
Power  Plant 
1.5MW*5 


Fig.  7.  On-shore  win  farms  in  Taiwan  ( Source :ITRI). 


a  b 

rum) 


Fig.  8.  The  linear  utility  function  (a)  Left  linear  membership  function  and  (b)  right  linear  membership  function. 
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spinning  turbine  blades).  Especially,  the  lower  frequencies  of  noise 
from  the  wind  turbine  are  judged  to  be  very  annoying  during  the 
night.  These  two  issues  are  directly  related  to  the  distance  from 
the  wind  turbines  to  residents  homes.  Thus  people  would  like  to 
maximize  the  distance  between  wind  turbine  installations  and 
their  dwelling  and  workplaces,  to  achieve  higher  satisfaction. 

But  due  to  the  costs  associated  with  construction  and  land 
acquisition,  and  considering  Taiwan's  dense  population,  it  is  often 
difficult  to  fully  meet  the  requirements  of  residents.  So,  finding  a 
balance  between  the  requirements  of  residents  and  construction 
considerations  is  the  most  important  issue.  For  sake  of  simplicity, 
three  basic  attributes  are  taken  into  consideration:  power  gener¬ 
ated  (MW),  investment  cost  (I),  and  social  acceptance  (S)  with  MF. 
Table  8  shows  the  data  for  MW  and  I  attributes.  Table  9  shows  five 
candidate  locations  for  the  installation  of  wind  turbines.  Table  10 
shows  the  social  acceptance  with  MF.  Based  on  (4.1)  and  Table  10, 


Table  8 

Alternatives  for  electric  generation. 


Alternative 

Variables 

MW 

I  (US$) 

Wind  turbine  1  (1  <  P  <  2  MW) 

1.98 

3,280,000 

Wind  turbine  2  (2  <  P  <  4  MW) 

*2 

3.92 

6,560,000 

Wind  turbine  3  (4  <  P  <  10  MW) 

*3 

9.12 

9,840,000 

Wind  turbine  4  (10  <  P  <  15  MW) 

*4 

14.55 

10,500,000 

Table  9 

Five  places  attributes. 


Places 

Distance  to  residents  house  (m) 

Land  cost  (US$) 

5 

300 

70,000 

6 

500 

90,000 

7 

600 

98,000 

8 

1200 

140,000 

9 

450 

80,000 

all  possible  levels  of  social  acceptance  (satisfaction)  for  wind 
turbine  locations  are  calculated,  as  shown  in  Table  11.  Four 
possible  wind  turbines  installations  (see  Table  8)  for  power 
generation  in  five  possible  locations  (see  Table  9)  must  be 
evaluated  in  order  to  satisfy  the  following  goals. 

-  (Gl)  Power  generated  must  be  higher  than  15  MW  and  the 
higher  the  better. 

-  (G2)  Investment  cost  must  be  limited  to  50  x  106$/year  and  the 
lower  the  better. 

-  (G3)  Social  acceptance  must  be  limited  to  the  level  of  50  and 
the  higher  the  better. 


Table  11 


Satisfaction  of  wind  turbine  installation. 

Alternative 

Place 

Value  of  MF  (satisfaction) 

Wind  turbine  1 

5 

0.125 

6 

0.375 

7 

0.5 

8 

1 

9 

0.3125 

Wind  turbine  2 

5 

0.052 

6 

0.263 

7 

0.368 

8 

1 

9 

0.21 

Wind  turbine  3 

5 

0 

6 

0.166 

7 

0.25 

8 

0.75 

9 

0.125 

Wind  turbine  4 

5 

0 

6 

0.09 

7 

0.18 

8 

0.727 

9 

0.045 

Table  10 

Wind  turbine  (noise  vs.  social  acceptance). 


Alternative  Distance  noise 


Wind  turbine  1  (1  <  P  <  2  MW)  [200  m,1000  m]  =  [40  db,30  db] 


Wind  turbine  2  (2  <  P  <  4  MW)  [250  m,1200  m]  =  [42  db,35  db] 


Wind  turbine  3  (4  <  P  <  10  MW)  [300  m,1500  m]  =[45  db,38  db] 


Wind  turbine  4  (10  <  P  <  15  MW)  [400  m,1500  m]  =[46  db,39  db] 


Social  acceptance 

[0,1] 


[0,1] 


[0,1] 


[0,1] 


Corresponding  MF 


(4.1)  right  linear  MF 


(4.1)  right  linear  MF 


(4.1 )  right  linear  MF 


(4.1)  right  linear  MF 


Graphs 


Satisfaction 
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Alternatives  Places 


The  achievement  model  without  adding  weights  to  each  goal  is 
formulated  using  the  proposed  model  MCGP-U  (see  the  Appendix 
for  details).  The  model  is  then  solved  using  LINGO  [28]  to  obtain 
the  solution  as  (MW,  I,  S)= (29.57,  34.26  x  10s,  0.4287).  Fig.  9 
shows  the  locations  for  each  wind  turbine  installation.  As  the 
result  the  following  goals  are  achieved:  social  acceptance  of 
42.87%,  power  generation  of  197%,  investment  cost  of  131.5% 
achieved.  Clearly,  the  problem  can  easily  be  solved  using  the 
proposed  MCGP-U.  In  order  to  increase  the  social  acceptance  rate, 
a  large  weight  is  added  to  G3.  Then,  the  problem  is  solved  again  to 
obtain  solution  as  (MW,  I,  S)=(16.53,  16.16  x  106,  0.6135).  Fig.  10 
shows  the  places  where  each  wind  turbine  should  be  located  for 
higher  social  acceptability.  As  seen  in  this  solution,  the  social 
acceptance  rate  increases  from  42.87%  to  61.35%  when  the  large 
weight  is  increased  to  G3.  This  shows  that  the  proposed  model 
provides  feasible  features  for  a  DM  to  deal  with  multiple  decision 
making  problems. 

5.  Conclusions 

Mitigation  of  climate  change  and  environmental  pollution  has 
increased  the  need  to  replace  fossil  fuel  energy  with  renewable 


energy  world-wide.  The  capacity  expansion  planning  problem  of 
renewable  energy  industry  is  an  important  issue  for  both  devel¬ 
oping  and  developed  counties  in  the  implementation  of  appro¬ 
priate  renewable  energy  policies.  To  deal  with  this  problem, 
Ramon  and  Cristobal  [9]  proposed  a  method  in  which  an  optimal 
outcome  of  plant  types,  numbers  of  plants  and  locations  can  be 
easily  determined.  In  this  study,  a  MCGP  model  is  proposed  to 
improve  the  usefulness  of  that  model  of  through  an  example  of 
locating  five  renewable  energy  plants  for  power  generation  in  five 
locations  in  the  autonomous  region  of  Cantabria,  in  the  north  of 
Spain.  Moreover,  the  normalization  method  is  introduced  to  the 
MCGP  for  helping  DMs  to  set  appropriate  weights  according  to 
their  preferences  or  policies.  As  a  result,  MCGP  can  be  used  as  an 
efficient,  effective  and  feasible  tool  for  DMs  when  dealing  with  the 
capacity  expansion  planning  problem  of  renewable  energy  indus¬ 
try.  In  addition,  social  acceptance  of  a  wind  turbine  installation  in 
Taiwan  is  also  considered. 


Appendix 

(MCGP-MSMS) 

Min  d1  +  d-2  +  e^~  ~^d2  T-  d4  +e^  ~\~d^  H-Cg-  T-  dg  + 

e6+d7 

s.t. — The  following  for  power  generated  goal  (the  case  of  the 


more  the  better): 

60, 959, 000Fi  8  + . . .  +  37, 770, 000  F7M  -  dt+  +  df  =yA,  (3.1) 

yA  -e+  +  ef  =  120  x  106,  (3.2) 

110xl06<y1<120  x  106,  the  sufficiently  large  value  for  power 
generated  is  120  x  106  W  (3.3) 

The  following  for  investment  (the  case  of  the  less  the  better): 
23, 425, 000Fi  8  + ...  +  47, 958, 000F7 14  —  +  d2  =y2,  (3.4) 

y2  -  e2  +  e2“  =80  x  106,  (3.5) 


80  x  106  <y2  <  100  x  106,the  sufficiently  small  value  for  investment  is  80  x  106 

(3.6) 

The  following  for  emission  avoided  goal  (the  case  of  the  more 


the  better): 

9, 649, 680F}  g  + . . . +4, 839, 548F7 14  —  +d2  =3^3,  (3.7) 

3/3  —  e^  =  19  x  106,  (3.8) 

18x106<3/3<19x  106,the  sufficiently  large  value  for  emission 

avoided  is  19  x  10  t  (3.9) 

The  following  for  jobs  created  goal  (the  case  of  the  more  the 
better): 

15Fi  g  + ... +  2OF7 14  —  d4h  +  d4  =  y4,  (310) 

yA-eX  +c4"  =90,  (3.11) 


70  <3/4  <  90,  the  sufficiently  large  value  for  jobs  created  is  90 

(3.12) 

The  following  for  operation  and  maintenance  costs  goal  (the 
case  of  the  less  the  better): 

37750Fi  g  + ...  +  75, OOOF7 14  —  d^  +d5“  =3/5,  (3.13) 

3/5  -e^~  +e$  =300000, 


(3.14) 
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300000  <y5  <  350000,  the  sufficiently  small  value  for  OM  costs  is  300000 

(3.15) 

The  following  for  total  distance  between  nodes  goal  (the  case  of 
the  more  the  better): 

6.5F8siiF9si  + ...  +  75F13siF145i  —  d6+  +  d6  =  y§,  (3.16) 

y6-e6++e6~  =4000,  (3.17) 

2910  <y6  <  4000,  the  sufficiently  large  value  for  total  distance  is  4000 

(3.18) 

The  following  for  social  acceptance  goal  (the  case  of  the  more 


the  better): 

7^1,8  +  •••  +  9F?  i4  —  d7  +d7  =50,  (3.19) 

f sosi  ~  F so\  —  FS02  —  FS02  —Fs04  —Fs05  —Fsoq  —Fso 7  =  0,  (3.20) 

F  8  si  +  f  9si  +  Ft  0  si  + f  1  lsi  +  Fusi  +  F 13  si  +  F 14  si  Fsiso  —  0,  (3.21) 

fl8  +^28  +^38  +^48  ~ Fssi  =  0,  (3.22) 

F29  +^49  +^69  —  Fgsi  =  0,  (3.23) 

^110 +^310 +^510 +^610  —  F-iosi  =  0,  (3.24) 


The  following  for  investment  (the  case  of  the  less  the  better): 

(3, 280, 000  +  700, 000)F15  +  ...  +  (10, 500, 000  +  800, 000)F49  -  dj  +  d2“  =  y2, 

(4.5) 

5, 000, 000  <  y2  <  1 0, 000, 000,  (4.6) 

y2-e2h  +e2  =  5,000,000,  (4.7) 

The  following  for  social  acceptance  (the  case  of  the  more  the 
better): 


0. 1 25F15  +  0.375F16  + . . .  +  0.045F49 
F 15  +F16  + ...  +F49 


d3  +d3  =y3. 


(4.8) 


0.5  <y3  <  1 


(4.9) 


y3-e3+  +e3  =1 


(4.10) 


The  following  for  one  wind  turbine  can  only  assign  one  place, 
and  vice  versa: 


X  Fij  <  1 ,  i  =  l,2,3,4 

j  =  5 


(4.11) 


ZFj<  1,  J  =  1,2,3, 4, 5. 


i  —  1 


(4.12) 


^111  "T  f  311  "Tf  511  “b  f  611  "Tf  711  —  fils/  =  0, 

(3.25) 

f  212  "b  f  512  "bf  712  f  12si  =  0, 

(3.26) 

f213  "b  f  413  +  f  613  f  13si  =  0, 

(3.27) 

fll4  +  f314+f714  —  fl4si  = 

(3.28) 

fsol  —  f  18  “ fllO  —  fll4  = 

(3.29) 

FSo2  ~  f  28  —  f  29  —  f  212  —  f  213  =  0, 

(3.30) 

fso3  —  f38  “  f310  “  f311  ~  f314  = 

(3.31) 

F so4  —  F  48  —  f  49  —  f  413  =  0, 

(3.32) 

FSo5  —  f510  “  f511  —  f512  —  0, 

(3.33) 

fso6  —  f69  —  f610  —  f611  —  f613  =  0> 

(3.34) 

fso7  —  f  711  —  f712  —  f714  =  0, 

(3.35) 

f  8  si  +  f  9si  +  f  1  Osi  +  f  1  lsi  +  f  12s/  +  f  13si  +  f  14si  =  5 , 

(3.36) 

where  Eqs.  (3.1)-(3.19)  correspond  to  the  goals  power,  investment, 
emissions,  jobs,  operation  and  maintenance,  and  social  accep¬ 
tance;  and  Eqs.  (3.20)-(3.36)  are  the  conservation  of  flow  con¬ 
straints  for  each  node. 

(MCGP-U) 

Min  dj  +  d^  +  6^  -bc^-  ~~b^2  ~~b^2  ~~b^2  ~^~^2 

s.t.  —The  following  for  power  generated  goal  (the  case  of  the  more 


the  better): 

1 .98Fi5  + 1 .98Fi5  + . . .  +  14.55F4g  —  d ^  +d|  =  jq  (4.1) 

10  <yt  < 30,  (4.2) 

y-i  ~  +ef  =30,  (4.3) 

The  sufficiently  large  value  for  power  generated  is  30  MW  (4.4) 


References 

[1]  Reserves,  Resources  and  Availability  of  Energy  Resources  2011;  German 
Mineral  Resource  Agency,  Hannover,  Germany. 

[2]  OPEC  share  of  world  oil  reserves  2011.  OPEC  2012. 

[3]  (http://www.ioga.com/Special/crudeoil_Hist.htm>. 

[4]  (http://www.infomine.com/investment/metal-prices/coal/5-year/). 

[5]  Hobbs  BF,  Meier  PM.  Multicriteria  methods  for  resource  planning:  an  experi¬ 
mental  comparison.  IEEE  Trans  Power  Syst  1994;9:1811-7. 

[6]  Afgan  NH,  Carvalho  MG.  Multi-criteria  assessment  of  new  and  renewable 
energy  power  plants.  Energy  2002;27:739-55. 

[7]  Terrados  J,  Almonacid  G,  Perez-Higueras  P.  Proposal  for  a  combined  metho¬ 
dology  for  renewable  energy  planning,  application  to  Spanish  region.  Renew 
and  Sust  Energy  Rev  2009;13:2022-30. 

[8]  Arnette  A,  Zobel  CW.  An  optimization  model  for  regional  renewable  energy 
development.  Renew  Sust  Energy  Rev  2012;16:4606-15. 

[9]  Ramon  J,  Cristobal  C.  A  goal  programming  model  for  the  optimal  mix  and 
location  of  renewable  energy  plants  in  the  north  of  Spain.  Renew  Sust  Energy 
Rev  2012;16:4461-4. 

[10]  Aras  H,  Erdogmus  S,  Koc  E.  Multi-criteria  selection  for  a  wind  observation 
station  location  using  analytic  hierarchy  process.  Renew  Energy  2004;  23: 
1383-92. 

[11]  Mladineo  N,  Margeta  J,  Brans  JP,  Mareschal  B.  Multicriteria  ranking  of 
alternative  locations  for  small  scale  hydro  plants.  Eur  J  Oper  Res 
1987;31:215-22. 

[12]  Beccali  M,  Cellura  M,  Mistretta  M.  Decision-making  in  energy  planning- 
application  of  ELECTRE  method  at  regional  level  for  diffusion  of  renewable 
energy  technology.  Renew  Energy  2003;28(12):2063-87. 

[13]  Cavallaro  F.  Fuzzy  TOPSIS  approach  for  assessing  thermal-energy  storage  in 
concentrated  solar  power  (CSP)  systems.  Appl  Energy  2010;87(2):496-503. 

[14]  San  Cristobal  JR.  Multi-criteria  decision-making  in  the  selection  of  a  renewable 
energy  project  in  Spain:  the  VIKOR  method.  Renew  Energy  2010;36:498-502. 

[15]  Goumas  M,  Lygerou  V.  An  extension  of  the  PROMETHEE  method  for  decision 
making  in  fuzzy  environment:  ranking  of  alternative  energy  exploitation 
projects.  Eur  J  Oper  Res  2000;123:606-13. 

[16]  Kaya  T,  Kahraman  C.  Multicriteria  renewable  energy  planning  using  an 
integrated  fuzzy  VIKOR  and  AHP  methodology:  the  case  of  Istanbul.  Energy 
2010;35(6):2515-27. 

[17]  Chang  C-T.  Multiple  choice  goal  programming  with  utility  functions.  Eur  J 
Oper  Res  2011;215:439-45. 

[18]  Chang  C-T.  Revised  multi-choice  goal  programming.  Appl  Math  Model 
2008;32:2587-95. 

[19]  Chang  Ching-Ter.  Multi-choice  goal  programming.  Omega  2007;35:389-96. 

[20]  Chang  C-T,  Ku  C-Y,  Ho  H-P,  Liao  C.  A  MCGP  decision  aid  for  homebuyers  to 
make  the  best  choice.  Qual  Quant  2011;45:969-83. 

[21]  Chang  Y-C,  Lee  N.  A  Multi-Objective  Goal  Programming  airport  selection 
model  for  low-cost  carriers'  networks.  Trans  Res  Part  E:  Logist  Transp  Rev 
2010;46(5):709-18. 

[22]  Liao  C-N,  Kao  H-P.  An  integrated  fuzzy  TOPSIS  and  MCGP  approach  to  supplier 
selection  in  supply  chain  management.  Expert  Syst  Appl  2011;38:10803-11. 


C.-T.  Chang  /  Renewable  and  Sustainable  Energy  Reviews  41  (2015)  379-389 


389 


[23]  Liao  C-N,  Kao  H-P.  Supplier  selection  model  using  Taguchi  loss  function, 
analytical  hierarchy  process  and  multi-choice  goal  programming.  Comp  Ind 
Eng  2010;58:571-7. 

[24]  Eyvindson  K,  Kangas  A,  Kurttila  M,  Hujala  T.  Using  preference  information  in 
developing  alternative  forest  plans.  Can  J  Forest  Res  2010;40:2398-410. 

[25]  Chen  Y-T,  Zheng  C,  Chang  C-T.  Efficiently  mapping  an  appropriate  thinning 
schedule  for  optimum  carbon  sequestration:  an  application  of  multi-segment 
goal  programming.  Forest  Ecol  Manag  2011;262:1168-73. 


[26]  Mahmoud  HB,  Ketata  R,  Romdhane  TB,  Ahmed  SB.  A  multiobjective- 
optimization  approach  for  a  piloted  quality-management  system:  a  compar¬ 
ison  of  two  approaches  for  a  case  study.  Comp  Ind  2011;62:460-6. 

[27]  Romero  C.  A  general  structure  of  achievement  function  for  a  goal  program¬ 
ming  model.  Eur  J  Oper  Res  2004;153:675-86. 

[28]  Extended  LINGO  Release  9.0,  LINDO  System,  Inc.  2004. 


